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ABSTRACT 

The water ice or snow line is one of the key properties of protoplanetary disks that determines the 
water content of terrestrial planets in the habitable zone. Its location is determined by the properties 
of the star, the mass accretion rate through the disk, and the size distribution of dust suspended in 
the disk. We calculate the snow line location from recent observations of mass accretion rates and as 
a function of stellar mass. By taking the observed dispersion in mass accretion rates as a measure 
of the dispersion in initial disk mass, we find that stars of a given mass will exhibit a range of snow 
line locations. At a given age and stellar mass, the observed dispersion in mass accretion rates of 
0.4 dex naturally leads to a dispersion in snow line locations of ^0.2 dex. Eor ISM-like dust sizes, 
the one-sigma snow line location among solar mass stars of the same age ranges from ^2 to ^5 an. 
Eor more realistic dust opacities that include larger grains, the snow line is located up to two times 
closer to the star. We use these locations and the outcome of N-body simulations to predict the 
amount of water delivered to terrestrial planets that formed in situ in the habitable zone. We find 
that the dispersion in snow line locations leads to a large range in water content. Eor ISM-like dust 
sizes, a significant fraction of habitable-zone terrestrial planets around sun-like stars remain dry, and 
no water is delivered to the habitable zones of low-mass M stars (< 0.5Mq) as in previous works. 
The closer-in snow line in disks with larger grains enables water delivery to the habitable zone for a 
significant fraction of M stars and all EGK stars. Considering their larger numbers and higher planet 
occurrence, M stars may host most of the water-rich terrestrial planets in the galaxy if these planets 
are able to hold on to their water in their subsequent evolution. 

Subject headings: planetary systems — protoplanet ary disks — planets and satellites: formation — 
planets and satellites: composition — stars: low-mass 


1. INTRODUCTION 

The leading explanation for the delivery of water to 
Earth is that water-bearing, asteroid-like bodies from be¬ 
yond the snow line were gravitationally scattered i nward 
and accreted by t he planet during its growth (e.g. Mor- 
bidelli et al.||2QQQ ). If the same mechanism also operates 
in extra-solar planetary systems, the water content of po¬ 
tentially habitable terrestrial planets that formed in situ 
depends on both the snow line location and the extent 
to which water-bearing materials were scattered inwards. 
Eor lower-mass M dwarf stars, N-body simulations indi¬ 
cated t hat terrestrial planets in the habitable zone would 
be dry (Raymond et al.|2007 Lissauer|2007), in contrast 
to planets that formed farther o ut in the disk and mi - 
grated to their curr ent locations (Ogi hara fc Ida||2QQ9). 
In a previous paper (Ciesla et aL|2UI5|), we explored how 


the inclusion of planetesimals and more comet-like bodies 
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in these simulations can enhance volatile delivery around 
these low-mass stars. In that work, the location of the 
snow line was a critical factor in determining how much 
water would be delivered to forming planets. Here, we 
revisit the location of the snow line based on recent pro- 
toplanetary disk observations, and explore its impact on 
water delivery to habitable zone terrestrial planets. 

In the Solar System, the location of the snow line at 
the time of planetesimal formation has been inferred from 
the transition b etween hydrous a nd anhydrous asteroids 
to be ^ 2.5 an (Abe et al. 2000). There is, however, a 
considerable uncertainty m this location if asteroi ds have 
been scattered into their current locations (e.g. Walsh 
et al. 2012 DeMeo fc Carry|[2014 ). There are no direct 


measurements of the location of a water snow line outside 
our solar system, hence we do not know if it is typical 
for sun-like stars. The Solar System snow line location 
of 2.5 an has been taken as a reference point in previous 
studies, and used to estimate the location around other 
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Figure 1. Location of the snow line for a range of stellar masses for a disk of small ISM-like grains (left panel) and larger grains more 
representative of protoplanetary disks (right). The green/red line, dark-shaded area and light-shaded area show the median, Tier, and ±2(7 
location of the snow line, respectively, calculated from the observed distribution of mass accretion rates using Eq. ^ 


stars of various masses (e.g. Raymond et al.||2007|). The 
exact location of the snow line will be deterniined by the 
combined ir radiative and viscous heating in the disk (e.g. 


Davis 


2005). If mass accretion occurs through the verti- 


cal extent of the disk, the mid plane temperatures at the 
relevant locations are set by the release of gravitational 
potential energy for mass accretion rates > 10“^^MQ/yr, 
and the location of the snow line is a strong function o f 
mass accretion rate (jGaraud & Lin|2007||0ka et al.|2011|). 

The dependence of mass accretion rate on stellar mass 
has been well-es tablished, following a roughly quadratic 


relationship (e.g.|Muzerolle et ah, 
Muzerolle et al. |2UU5| Natta et al.| |2 


Calvet et al. 2004 
but these re¬ 
lations have not been used to calculate the location of 
the snow line around low-mass stars. On top of that, 
there is an intrinsic dispersion in mass accretion rate s of 
^ 0.4 dex around this relation (Alcala et al. 2014). If 
this intrinsic dispersions reflects a range m i nitial disk 


masses, wh ich is of the same order, 0.5 dex (Armitage 
erni[2M| , we expect an associated dispersion m snow 


line locations (Figure [^, leading to enhanced or reduced 
water delivery among stars of similar mass. 

In addition to the mass accretion rate, the location 
of the snow line in a protoplanetary disk aro und a sun- 


like star is very sensit ive to the dust opacity (Min et al. 
2011||Okaet al.|2011 ). While previous work on the snow 


line location in the solar ne bula had ma i nly used small 


sub-m icron sized grains (e.g. Davis 2005 Garaud & Lin 


2007), there is abundant observational and theoretical 
evidence that grains in protoplanetary disks quickly grow 
to millimeter and centimeter sizes (e.g. Testi et al.|2014), 
depleting the amount of sma ll grain by orders of magni- 
tude (e.g. jFurlan et al.||2006|). Lower dust opacities trap 
viscous heat less efticiently, resulting in a cooler disk mid 
plane with a snow line closer to the star, enhancing water 
delivery. 

In this work, we explore the range of snow line loca¬ 
tions that may be expected for stars of different masses, 
factoring in variations in disk mass accretion rate and 
dust opacities that have been inferred for real disks (^ 
and how this would impact the water content of terres¬ 
trial planets in the Habitable zone (^. 


2. SNOW LINE LOCATION IN A VISCOUS DISK 

The location of the snow line, Rsl, in a steady-state 
viscous protoplanetary disk, assu ming it is optically thick 
to its own radiation, is given by Min et al. (2011), their 
Eq. (11): 


( 128 w'^kbasBTicefa J 


( 1 ) 


where /i is the mean molecular weight, rup the proton 
mass, G the gravitational constant, M the mass accre¬ 
tion rate, the stellar mass, kr the Rosseland mean 
opacity, the Boltzmann constant, ctsb the Stefan- 
Boltzmann constant, Tice the temperature where water 
ice condenses, / the gas-to-dust ratio, and a the tur¬ 
bulent mixing strength. This equation is based on the 
est imate of t he mi d plane temperatur e of a viscous disk 
by Hubeny (1990). Min et al. (2011) showed that this 
equation is in good agreement with the snow line com¬ 
puted through detailed radiative transfer modeling for 
a sun-like star. Even though the role of irradiation be¬ 
comes larger when using the pre-main sequence luminos¬ 
ity and when considering lower mass disks around lower 
mass stars, we have verified in Appendix |A] that viscous 
heating indeed dominates the thermal bud^t of the disk 
mid plane for the range of stellar masses and disk mass 
accretion rates explored in this paper. 

The ice sublimation tem perature is thought to li e be¬ 
tween 150 and 170 K (e.g. Podolak fc Zucker]|2004 ), de¬ 
pending on the water vapo r pressure in the disk. The 

(|2011|) include a 


Min et al. 


radiative transfer models in 
detailed treatment of solid sublimation as described in 


(Kama et al. 2009). Min et al. (2011) show that, for 


a large range of mass accretion rates — and hence sur¬ 
face densities and partial vapor pressures — the location 
of the snow line is well predicted by assuming a single 
sublimation temperature of Tice = 160LC. We therefore 


do not solve the radial diffusion equation (e.g. Ciesla & 


Cuzzi 2006), as we expect the variations of the partial 


water pressure by inward drift of evaporating ices to be 
within the tested range and not lead to large variations 
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tet^ 

log M 

^SL D-0-) 


^SL 



small grains 

large grains 

MS 

pre-MS 

[Mq\ 

[logMo/yr] 

[au] 

[au] 

[au] 

[au] 

1.0 

-7.5 ±0.4 

3.3 (2.2...5.0) 

1.6 (1.1...2.4) 

2.5 

1.35 

0.8 

-7.7 ±0.4 

2.6 4-7...3.9) 

1.2 (0.82...1.9) 

1.23 

1.15 

0.6 

-7.9 ±0.4 

1.9 4-2...2.8) 

0.9 (0.59...1.3) 

0.61 

0.95 

0.4 

-8.3 ±0.4 

1.2 5-78...1.8) 

0.6 (0.37...0.84) 

0.32 

0.8 

0.2 

-8.8 ±0.4 

0.5 (0.36...0.81) 

0.2 (0.17...0.38) 

0.16 

0.5 


Table 1 


Mass accretion rates and location of the viscous snow line as a function of stellar mass used in this work, including one-sigma ranges. For 
comparison, the last two column shows the s now line locations calc ulated from t he stellar luminosi ty used in previous work: MS 
(main-sequence, [Raymond et al.|2007| ) and pre-MS ( [Ciesla et al.|2Ql^ . 


in the location of the snow line. 

Grouping all constants together, taking fi = 2.3, we 
obtain the more manageable expression 


i?SL = 2.1 au ( ^ ) 


M 


MeJ \ 10-8M©/yr 

2/9 / n \ -2/9 


I^R 


770 cm^/g 

-2/9 


f—)■ 

Vo.oi/ 


' ( 2 ) 

100 J 

-10/9 


160 K 


In the remainder of the section we explain the origin and 
applicability of the parameters in equation that lead to 
the range of snow lines shown in Figure 

2 . 1 . Disk mass accretion rate 
The mass accret ion rate scales roughly with stellar 


2004 


squared (e.g. Muzerolle et al. 2003[ 

Cal vet et al. 

Muzerolle et al. 2UU5| hJatta et al. 

| 2 UU 6 |), with 


an observed dispersion of up to an order of magnitude 
around this relation. We will use the observed mass ac¬ 


cretion rates from Alcala et al. (2014) since these have 
the sm alles t observed dispersion, and we will explain in 
section 2.2 how we treat this dispersion. 

The median mass accretion rate for the ^ 3 ± 1 Myr 


old (Comeron 2008) Lupus star forming region, over a 
ra nge of stellar masses from 0.03 Mq to I.OMq, is given 
by Alcala et al. (2014): 


-^Lupus(-^*) = 5.6 • 10 


-9 


M* 

Mcr 


-I 1.81 


Me/yr. (3) 


We derive a generic expression for the median mass 
accretion rate as a function of time by appending this 
equation with a factor to take into account the age 
of the cluster: 


1 - 3/2 


t) = 3 • 10 ^ —— M 0 /yr. 

^ ^ _M©J [MyrJ 

(4) 

This accretion rate is consistent with the derived ac¬ 
cretion rate for the 1 Myr ol d Taurus region of ^ 
lO“^M 0 /yr for ^ 0.5 M© stars ( Hartmann et al.flQQS ), 
and matches well with the observed dec ay of the mass 
accretion ra te between 0.3 and 30 Myr (Sicilia-Aguilar 
et al. 2010 their Figure 2 ). The exponent of 77 = 1.5 
comes irom a viscously evolvin g disk with fixed temper¬ 
ature ([Hartmann et al. ||2006|), and matches well with 
the analytical model of 


Chambers (2009) that includes 


a changing disk temperature {g = 20/13 = 1.54) which 
is more similar to the steady sta te accretion disk m odels 
employed in deriving Eq. fl] by [Min et al.j (|2011|) . Eq. 

is only valid at time scales much longer than tne vis¬ 
cous time scale, and starts significantly over-predicting 
the mass accretion rate at t < 0.3 Myr for a = 0.01 


compared to the Chambers (2009) disk model. 


2 . 2 . Range of snow line locations 

A dispersion in disk mass accretion rate at a given age 
will lead to a dispersion in snow line locations. Observed 
mass accretion rates show a large dispersion at any given 
stellar mass. Part of this dispersion may be att ributed to 


variations i n the accretion flow onto the star (Costigan 
and use of secondary tracers of accretion 


2012 ), 


et al. 

such as emission J 
dispersion may reflect an age spread and variati ons in 


. lines ( jRigliaco et al.| 2 QT 2 ). An intrinsic 


Hartmann et al. 


the initial disk mass (e.g. 
tightest con straint on the intrin sic dispersion currently 


1998). The 


comes from Alcala et al. (2014). Using the X-shooter 
spectrograph, the authors compute the mass accretion 
rates using the UV excess emission, a direct tracer of ac¬ 
cretion, in the Lupus star forming region over a range 
of stellar masses from 0.03 M© to l.OM©. The disper¬ 
sion in mass accretion rates at a given stellar mass is 
characterized by a log-normal distribution with a stan¬ 
dard deviation = 0.4 dex, a factor of two lower than 
in previous works. We assume this dispersio n reflects a 
range of ini tial disk masses as proposed by jHartmaim] 


^tal ([1998 ). 

I'his assumption is supported by two results. Eirst, the 
diversity in protoplanetary disk pr operties, in particular 


the outer radius-disk mass relation ( Andrews et al.| 20 TQ ) 
and the exponential decay of the disk fraction versus time 
(e.g. Mamajek||2009) argue for an intrinsic spread in ini¬ 
tial disk massed Since the inferred spread in disk masses 
from the time-de pendent disks fractio n is of the same 
order, ^ 0.5 dex ([Armitage et al.||2003|), we take the ob¬ 
served dispersion m mass accretion rate to be the intrin¬ 
sic one. Second, although accretion rates measured from 
spectral lines are known to be variable on time scales of 
days to years, this variation is smaller than the observed 
dispersion and driven by the r otational modulation of the 
accretion flow onto the star (Costigan et al. 2014), and 
does not reflect an intrinsic variation in the disk mass 
accretion rate. It is possible that disks show larger vari¬ 
ations at time scales of decades or longer but smaller in 
magnitude than EUor outbursts, or that it reflect a large 
age spread in the Lupus cloud, in which case the intrinsic 
scatter could be smaller. 

With the assumption that the observed spread in mass 
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accretion rates is the intrinsic one, the probability distri¬ 
bution of the snow line R around the median, i?sL (Ed- 
1^, is given by a log-normal distribution: 


/sl(-R,-Rsl) = ^ „ -exp 

V27r R an 


(log-R-log-RsL)^ 

24 


( 5 ) 

with a standard deviation cjr = 4/9 * ~ 0.18 dex. 

This distribution is shown as the shaded regions in Fig. 

m 


2.3. Dust opacity 
Since a major source of uncertainty in the calculated 


locat ion of the snow line is the dust opacity (Min et al. 
2011), we explore two limiting cases: the first is where 


the grains are ISM-like in size as in most previous work, 
with a Rosseland mean opacity at temperatures relevant 
for water i c e con densation of kr = 570 cm^/^, as in 
Min et al. (2011); the second is where dust growth is 
assumed to occur, resulting in a size distribution from 
micron to centimeter-sized grains, yielding an opacity of 
Hr = 20 /g. Such a lower opacity will also increase 

the amount of solids availab le for planetesima l formation 
interior to the snow line, see|Min et al. (2011) for details. 
Because the dependence of the snow line location on this 

opacity is relatively weak (Rsl Eq[^, the snow 

line moves in by about a factor of two between the two 
cases considered here (left and right panels of Fig. [^. 
Table lists the mass accretion and corresponding snow 
line locations for the range of stellar masses for the two 
opacities considered here. Figure shows the median 
snow lines for the same stellar mass range , compared 
to the recen t habitable zone estimate from |Kopparapu| 
et al. (2013). Throughout this paper, we use a gas-to- 
dust ratio / = 66 based o n a condensation sequence at 
solar metallicity, following Min et al. (2011). 


2.4. Snow line location 

For a solar mass star with small grains, the one-sigma 
range of snow line locations encompasses that of the So¬ 
lar System, from ^2 to ^5 au . Its m edian location is 
consistent with the Chambers (2009) disk model for a 


disk of the same mass accretion rate, age, and opacity. 
Due to the relatively weak scaling with mass accretion 
rate the two-sigma range of snow line location 

changes by a factor five, even if M changes by a factor 
of 40. For the large grain case, the one-sigma locations 
range from 1.1 to 2.4 au, roughly encompassing the hab¬ 
itable zone. In the context of these models, if the solar 
system snow line was located at 2.5 au at 1 Myr, it is 
a one-sigma out-lier, and the majority of exoplanetary 
systems may form with a snow line closer to the star. 

Around lower mass stars, the snow line moves closer to 
the star as Rsl oc (Eqs[^[^. The snow line around 

a star is closer in compared to that inferred by 

Kennedy & Kenyon (2008), who use a linear scaling of the 
mass accretion rate with stellar mass which is no longer 
supported by observations, leading to a hotter disk with 
a more distant snow line. Compared to the habitable 
zone, that scales with stellar mass roughly as Rhz 

MS oc relative gap between snow line 

and habitable zone widens (See also Lissauer|2007 ). For 



Figure 2. The location of the snow line in au around a solar mass 
star, for a range of the turbulent mixing strength a and the time t 
at which the snow line is “frozen in” the planetesimal population. 
Our choice of a = 0.01 and t = 1 Myr is shown by the central cross. 
In general, lower turbulent mixing strength and earlier times move 
the snow line outward. Towards higher turbulent mixing strength 
and later times it moves inward. 

the small grain case, the gap between the median snow 
line and the habitable zone is about a factor of five. For 
the large grain case, this is only a factor of two, and the 
two-sigma tail of the snow line distribution overlaps with 
the habitable zone. 


2.5. Disk viscosity and time of planetesimal formation 

Throughout this paper, we assume t = 1 Myr and 
a = 0.01 for setting the location of the snow line. Both 
quantities are not well constrained observationally. The 
strength of viscosity in protoplanetary disks is inferred 
to be of orde r a ^ 0.01 based on disk masses, sizes, 
and lifetimes (Hartmann et al.||1998). The time that de¬ 
fines the location of the snow line depends on both the 
planet formation time scale and the extent to which (icy) 
grains continue to be incorporated into larger bodies as 
the snow line moves in. The theoretical planetesimal for¬ 
mat ion__yin^_scaJT_a^^^lA^is of order 0.01 — 0.1 
(e.g. Wetherill & Stewart|1993), and t he earliest differen ¬ 
tiated bodies formed within < IMyr (jKleine et al.|2009|), 
suggesting an early epoch of planetesimal formation. On 
the other hand, the age of the chondrules in chondritic 
parent bodies indicate planetesimal form ation continued 
to tak e place over a 3 Myr time scale (e.g. | Johansen et al. 
2014|). If water vapor remains present in the disk, an in¬ 
ward moving snow line may give rise to rapid, bu t late 
icy planetesimal formation (Kennedy et al. 2006|). The 
age of t = 1 Myr is a compromise between two extremes. 

To give the reader an idea of how these uncertainties 
impact the location of the snow line. Figure shows how 
the location of the snow line varies around a sun-like star 
with a median mass accretion rate following equation 
^ Earlier times and lower turbulent mixing strength 
correspond to higher surface densities and mass accretion 
rates that lead to a hotter mid plane and a more distant 
snow line. Higher mass accretion rates and later times 
have the opposite effect, leading to a colder, optically 
thinner disk with a closer-in snow line. In this regime, 
irradiation may become important as a heating source in 
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Figure 3. Planetary water abundances for a 0.6Mq star, com¬ 
puted for three different snow line l ocations (median a nd ±2(7). 
Compare with CC case in figure 6 in |Ciesla et al.| ( [2Q15[ ). 

the mid plane. 

3. IMPLICATIONS FOR WATER DELIVERY 

To explore how the range in snow line locations may 
impact the water content of terrestrial planets in the 
habit able zone, we use t he set of N-body simulations 


fro m iCiesla et ah (|2Q15[ ) . This approach follows that 


of lU'Bri e n et al.T ( 
Lissauer (|2UU7|), f< 


Raymond et ah (2007) and 


ti007| , I or forming terrestrial planets via gi 


ant impacts as in the solar system. The simulations are 
started from a disk with equal mass planetary em¬ 
bryos (1/20 Earth mass) and planetesimals (1/20 embryo 
mass), spaced such t hat they define a 1/R power law in 
surface density, as in|0’Brien et al.| (|2006|), between 0.5 
and 4.0 au. We do not include giant planets, but note 


that they m ay be rare around low-mass stars (Johnson 
et al. 2010|) and water delivery can take place in their 


absence dQuintana & Lissauer |2014|). Four simulations 
are run with near-identical initial conditions to take into 
account stochastics. These are repeated for lower stellar 
masses (0.2, 0.4, 0.6 and 0.8 Mq) by scaling all masses 
(total, embryo, and planetesimal) to the stellar mass and 
adjusting the inner and outer edge to encompass both 
snow line and habitable zone. After 200 Myr, a number 
of terrestrial plane ts has formed in the inner regions, see 
Ciesla et al. (2015) for details. 

T'he water content of these planets depends on the 
composition of accreted planetary embryos and planetes¬ 
imals, which in turn depends on their starting locations. 
We assume that the division between icy and rocky build¬ 
ing blocks is located at the snow line, i.e. that plan¬ 
etesimals and embryos do not move radially during their 
growth process, and that the transition from pebbles to 
planetesimals and embryos is short compared to the time 
scale at wh ich the snow line moves inward (e.g. Carreraj 
et al.||2015|) . Assuming the water content of icy building 
blocks is low (5%, similar to carbonaceous chondrites), 
the extra mass outside of the snow line does not affect 
the dynamics of the system, allowing the snow line to 
be inserted a posteriori. By keeping track of the start¬ 
ing locations of planetary building blocks that comprise 
the final planets, their water content can be calculated 
for different snow line locations from a single simulation. 


Figure gives an example of this approach, which shows 
the water content of planets in the 0.6 Mq simulations 
calculated from three different locations of the snow line. 

By calculating the water content of planets in the hab¬ 
itable zone for a range of snow line locations from the 
simulations, and assigning to each snow line a probabil¬ 
ity based on the observed dispersion in mass accretion 
rates (Eq. see also Fig. [^, we calculate the proba¬ 
bility distribution of the water fraction of planets in the 
habitable zone. These are shown as cumulative proba¬ 
bility distributions in Fig. 

We take a water mass fraction of 0.1% as a reference 
point fo r early Earth, b etween 0.03-0.1% currently in the 
mantle (Lecuyer 2013|) and as much as 1% in the past 


according to Abe et al. (|2000|). In the simulations with 
sun-like stars, for a snow line set by small grains, 20% 
to 40% of planets in the habitable zone remain dry, as a 
fraction of stars has a mass accretion rate high enough 
to put the snow line too far away from the habitable 
zone ( 4-5 au) to allow water delivery via giant impacts. 
The fraction of dry planets is higher in simulations with 
lower-mass stars, reaching 40% to 60% for 0.6 and 0.8 
solar mass stars, and close to 100 % for < 0.4 solar mass 
stars, in line with previous results (Raymond et al.|2007 
Ciesla et al.||2015|). The one-sigma range of water trac- 


tions spans more than an order of magnitude, signifi¬ 
cantly larger than the stochastic dispersion in the N-body 
simulations (a factor of ^2 to ^5). 

Using larger grains for calculating the snow line loca¬ 
tion significantly changes the results. Because the snow 
line is located closer to the star, water delivery is more 
efficient: the fraction of dry habitable-zone planets in our 
simulations decreases to zero for 0.8 and 1.0 solar mass 
stars, and to 10% for 0.6 Mq. A small fraction of hab¬ 
itable zone planets around low-mass stars does receive 
earthlike amounts of water, around 20% for 0.2 Mq and 
0.4 Mq. 

4. DISCUSSION AND CONCLUSIONS 

We derived the fraction of water-bearing terrestrial 
planets in the habitable zone by assuming they are 
formed in-situ. There is, however, growing observa¬ 
tional evidence that migration of planets or their build¬ 
ing blocks plays an important role in the formation of 
super-Earths and mini-Neptunes at short orbital periods 


(iHansen fc Murr^|2012| [Swift et ahll^ 
Cossou pTlt l llVIulders et al. |:^ul5X Wh 


2013 Raymond & 


|2U15|. Whether this mecha¬ 
nism is' able to form smaller, earth-sized planets farther 
out remains an open question. If this were the case, the 
population of terrestrial planets in the habitable zone 
might be a mixture of planets formed in-situ and through 
migration. As mi grating planets typically form beyond 
the snow line (e.g. Cossou et al.|2Q14 Izidoro et al.|2Q14 ), 
they are expected to have water masVfract ions of order 
50% (jOgihara & Ida|2009t |Tian & Ida|2015|. Hence, the 
fractions of wet terrestrial planets quoted in this paper 
can be considered as lower limits. 

An additional concern for the habitability of habitable 
zone planet around lower mass stars is that water loss 
may may be more efficient as t hese planets fo r m hotter 
than those around sun like stars (Lissauer 2007|) . Ramirez 
& Kaltenegger (2014) show that duing the bright pre- 
mam sequence phase of low-mass stars, stellar fluxes are 
high enough to trigger a runaway green house effect that 
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Water Distribution in Habitable Zone 



Water Distribution in Habitable Zone 



Figure 4. Fraction of planets in the habitable zone with at least a given mass fraction of water, for stars of different masses. The left and 
right panel use a snow line calculated with ISM-like and larger grains, respectively. 


leads to enhanced water loss. Tian & Ida (2015) come 
to a similar conclusion by modeling water los tnrough 
hydrodynamc escape in a planet population synthesis 
model. For a planet to remain habitable, water has to 
arrive later during the pre-main-sequence evolution of 
the star. The timing and delivery mechanism of water 
may be crucial here: In our simulations, water delivery 
persists over ^ 100 Myr time scales. In addition, water 
is not delivered in the form of ices that would be di¬ 
rectly deposited in the atmosphere. Rather, it arrives as 
hydrated minerals that would enter the atmosphere later 
via vulcanic outgassing. A detailed study of arrival times 
and atmospheric release of water may be neccessary to 
address whether water can be deliverd late enough and 
in sufficient quantities to avoid escape. 

Despite water delivery to low-mass planets in the hab- 
i table zone being less efficient than aro und sunlike stars 
(Lissauer 2007 Raymond et al. 2007), we show that a 
small fraction niay still receive earthlike amounts of wa¬ 
ter due to the dispersion in snow line locations. When 
taking into account the larger number of M stars with 
respect to G stars (factor ten), their higher p lanet- 
occurrence rates (factor two, Mulders et al. 2015|), and 
i ncreased water de livery from more comet-like icy bodies 
(Ciesla et al. the majority of water-rich terrestrial 

planets may still be found around low-mass M stars. 
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Figure 5. Ratio of the snow line location calculated from radiative transfer models with respect to the analytical approximation (Eq. [^. 
The blue shaded region indicates a 10% error. The analytical approximation shows good agreement with the radiative transfer models over 
the ±4 (7 range of mass accretion rates for stellar masses > 0.6Mq. At lower stellar masses, irradiation leads to a larger snow line at mass 
accretion rates 3 and 2 sigma below the median mass accretion rate for a stellar mass of 0.4Mq and 0.2Mq, respectively. 
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APPENDIX 

A. THE ROLE OF IRRADIATION IN DETERMINING THE SNOW LINE LOCATION. 


Throughout this paper^e have assumed that the midplane temperature at the location of th e snow line is determined 
by viscous heating (Eq. fT|, and irradiation is negligable. For sunlike stars, Min et al.|(2 Qll|) have shown that this is 
a goo d approximation for the mass accretion rates considered in thi s paper (bee also [Garaud fc Lin||2QQ7 Oka et ah 


2011), assuming accretion is not layered (e.g. Lesniak fc Desch|[2Qll ). Around lower mass stars, lower mass accretion 


rates and a brighter pre-main sequence phase may lead to a larger role of irradiation. To verify the assumption t hat 
equation [I] also holds for low-mass stars, we compute the original radiative tra nsfer model from|Min et al.|(2011) for 
lower-mass stars. For the stellar photosphere, we use the Baraffe et al.| (1998) evolutionary tracks at t = 1 Myr. We 
explore a ±4cr range in mass accretion rates for the stellar masses used in this paper. Figure!^ shows the location of 
the snow line (i.e., the location where half the water is condensed into ice) divided by the predicted location from Eq. 

m 


For stellar masses larger than O.4M0, the location of the snow line is predicted accurately down to two sigma below 
the median within the precision of the radiative transfer model 5%). At smaller mass accretion rates, irradiation 
becomes imporant and equation [^starts underpredicting the location of the snow line. For the 0.2 Mq case, irradiation 
has a larger influence, and the predicted location of the snow line starts deviating at two sigma below the median. 






















